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We study the topological properties of the quantum phase space current in the 
Husimi representation, focusing on the dynamical differences, induced by these prop¬ 
erties, between the quantum and the classical flows. We show that the zeros of the 
Husimi function are stagnation points of the current and have a nonzero topologi¬ 
cal charge. Due to overall charge conservation, new stagnation points with opposite 
charge appear in pairs in the Husimi current and they have important roles in dy¬ 
namical processes. As an example we show the topological effect of the zeros in the 
transmission rate of particle through a potential barrier. 



2 


I. INTRODUCTION 

Motivation. Since the early days of the quantum theory physicists have looked for ways 
to represent quantum states as probability distributions in the phase space, a procedure 
in which information about coordinate and momentum representations is encoded in one 
single function M- This type of representation, thoroughly employed in classical statis¬ 
tical mechanics, is nowadays widely used to study quantum systems with applications in 
experimental photonics, quantum information, numerical semiclassical methods and other 
applied and theoretical fields. 

In general, functions Q representing probability densities in the phase space satisfy the 
continuity equation 

^ + divJ = a, (1) 

where the divergence is taken with respect to the coordinates x and momenta p. This 
equation defines the probability current J and the generation rates cr for Q, which can be 
either positive or negative. In classical mechanics cr = 0 due to probability conservation, 
and the equations of motion for the coordinates determine the current |6j. For example, in 
an one-dimensional system 



where the dots indicate total derivative with respect to time. Classical solutions of the 
equation ([TJ) can be obtained from an initial function Q and are given in terms of the 
trajectories governed by the equations of motion: Q(x,p]t ) = Q (xo,Poi 0); where Xq and 
Pq are the initial conditions of the trajectory ending at x and p in the time t. These 
classical probability functions have two properties: their marginal distributions are the 
physical coordinate and momentum distributions and they are positive semi-definite, Q > 0. 

The quantum mechanical analog of this formulation uses phase space functions to repre¬ 
sent states but the general solutions have different properties from that of the classical ones. 
Firstly, it has long been known that any quantum function violates at least one of the two 
classical properties stated before. For example, the Wigner function lam has as marginal 
distributions the wavefunctions of the state in the coordinate and momentum representa¬ 
tions but it is not positive, and the Husimi function M is positive, but it does not have 
the wavefunctions as its marginal distributions. Secondly, quantum solutions are not guided 
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by trajectories because it is impossible to assess both coordinate and momentum for a given 
test particle, even though the quantum functions satisfy the continuity equation and have 
well a defined probability current mm- Since authentic quantum trajectories do not exist 
as in the classical dynamics, in order to characterize both dynamical regimes in the phase 
space on the same grounds we need to attain the analysis on the probability flow itself. 

Quantum phase space flow. Continuity equations in quantum mechanics cannot be de¬ 
rived from equations of motion. Instead they are built by casting the von Neumann equation 
for the dynamics of the states in some representation, in which the density operator of the 
state is mapped to the function Q in the phase space, where a suitable current J is obtained. 
In recent years attention has been drawn to the topological structures the quantum current 
generates and the relations to its classical counterpart. So far it has been shown that the 
Wigner current has a strict topological order in the dynamics of its stagnation points ins. 
and in a recent work the same kind of order was explored for the Husimi function m 

The Husimi function is obtained by averaging the density operator in a basis of Gaussian 
coherent states and it can be interpreted as the probability of measuring the position and 
momentum of the quantum system within an uncertainty area centered in the basis state. 
This function inherits the analiticity of the coherent states, which restricts the dynamics of 
the stagnation points. In a previous work HH we presented a detailed demonstration of the 
construction of the Husimi current for one-dimensional systems, and we showed that the 
current adds new topological structures to the dynamics when compared with its classical 
counterpart. We also conjectured that these new structures emerge because every zero of 
the Husimi function is also a saddle point of the flow. In this paper we prove this conjecture 
and we explore the effects of this topological structures in a toy model. 

Specifically, we consider transmission through a gaussian barrier and show the connection 
between the classical and quantum transmission coefficients and the position of the zeros 
relative to the peaks of Husimi function and the classical energy levels. It has already 
been observed that the emergence of zeros in the phase space is a signature of quantum 
interference m fl3j : our work adds another layer of comprehension of this signature by 
analysing a dynamical feature of the zeros. 

Structure of the paper. In section 2 we present a summary of our previous work m- 
Section 3 contains the demonstration of the conjecture for the two-dimensional phase space. 
In section 4 we use the Gaussian barrier as a toy model to analyse the dynamical effect of 
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the topological structures. Section 5 contains the final remarks. 

II. THE HUSIMI FLOW 

The coherent states | z) of the harmonic oscillator of mass m and frequency c o are defined 
as the eigenstates of the annihilation operator, 

a |*) = z\z), 

where z is a complex number that labels the eigenvalue. The normalized coherent states 
can be written as 

| z) = e.- 2 " z/ V st |0), 

where z is the complex conjugate of z and |0) is the ground state of the harmonic oscillator. 

Throughout this paper we use a bar to denote the complex conjugate. The annihilation and 

creation operators, a and cd, are given by 

^ x . p x . p 

a = -b i —, a—- i —, 

where a x = yjh/2muj and a p = \Jhmoj /2 are the coordinate and momentum widths of the 

ground state, respectively 0, II- The coherent states are minimum uncertainty localized 

Gaussian wavepackets centered at x — (z|a;|z) = a x D\c (z) and p = (z\p\z) = a p 3m(z). In 

terms of x and p we can rewrite 0 as 

x p x p 

z = -b i —, z — - i — 

@ X 

The change of variables (x,p) i—> ( z,ihz ) is a canonical transformation and due to the 
scaling by cr x and a p the dimensionless z coordinate have characteristic size of 1/which 
is important when expansions in powers of H are needed. 

The identity operator in the coherent state representation is 

1 = J A 2 z\z)(z\, 

where the d 2 z = d^d^/27rz = dxdp/2-Kh is the displacement invariant volume of the phase 
space. A quantum state represented by the density operator 'p can be mapped to a phase 
space function by means of its average over coherent states: 



Q(z,z) — tr (p\z){z\). 


( 3 ) 
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This is the Husimi function of the quantum state, also called its Q-symbol. If the state is 
pure the Husimi function is the squared modulus of the wavefunction and is non-negative. It 
represents the minimal uncertainty probability of measuring the position and momentum of 
the state. From the definition of |z), the wavefunction in the coherent state representation 
can be written as 

(z, z) = {z\if>) = e~ gz/2 9 (z ), (4) 


where 9 (z) = (0\e za \ip) is an analytic function of z Pi- Similarly, ^{z,z) = z) = 
e~ zz / 2 9 (z). This factorization of the wave function will be important in the next section. 
The dynamics of the density operator is governed by the von Neumann equation, 




( 5 ) 


where H is the Hamiltonian operator and [•, •] is the commutator. In order to represent the 
dynamics in the phase space, we define a Hamiltonian function for H, 


H(z,z) = tr (h\z){z\) , (6) 

which is an average of the operator over the coherent states. If the Hamiltonian operator 
can be written as a normal ordered power series in a and a', namely H = Yhm h mn cd m a n , 
then its averaged function is H = h mn z m z n . With this definition, equation (|5j) can be 
rewritten as a differential equation in the phase space for the Husimi function: 

. r dQ x- _ m fd y 
ih— = 2^h mn z + Q-c.c., 

m,n ' ' 

where c.c. stands for the complex conjugate of the expression immediately preceding it. 
Taking the limit h —y 0 the right hand side of this equation reduces to the Poisson bracket, 
and the classical dynamics is obtained. In this limit the evolution of a phase space function 
can be further written as a continuity equation: 

(JO dJ c i dJ c i 

— H-- H-- = 0 

dt dz dz 

where J c i (z,z) = is the classical probability current. This classical dynamics is based 

on the existence of trajectories, guided by the equation of motion z — that carry the 

Q function over the phase space. In a similar fashion we can transform the full quantum 
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dynamical equation into a sourceless continuity equation if the Hamiltonian operator is 
Hermitian: 


dQ dJ dJ 
~dt + ~dz + ai = ’ 

where J (z, z ) is the quantum probability current and is given by 


(7) 


1 ^ mi ^ } ^ h mn (—l) fc m\n\ d 1 - 1 


J = V _ /_ m -k-ln-kft\ 

m,n =0 fc=0 


ih 


Z-j' fc!Z! (m — k — l)\ (n — k)\ dz l ~ 


( 8 ) 


The lowest order h term on the right hand side of (|8j) is obtained with l = 1 and k — 0, 
and this term retrieves the classical current J c i. In this way, the quantum current can be 
separated into the classical current plus higher order corrections in h. If the Hamiltonian is 
not Hermitian, h mn h nm , a source term 

a = (H-H) 

in v ' 


is added to the right hand side of ([?]). This term accounts for the absence of norm conser¬ 
vation of the Husimi function. In this paper we will only consider Hermitian Hamiltonians. 

Simplification of J. The expression (J8| can be further simplified by changing the sum¬ 
mation indexes in the following way: 


oo oo min(ra,n) m—k oo oo oo oo 

EE E E-EE E E- 

m= 0 n=0 fc=0 1= 1 fc=0 Z=1 m=k+l+ 1 n=k 

Expanding the derivatives in z contained in (J8]) and gathering conveniently the terms we 
obtain 

1 y. tf-'Q y. (-1)* d^H 

ih dz l ~ l (k -g IV dz k+l dz k ’ 

1=1 k =o v ’ 

which is a more compact and manageable expression than the original ([8]). Our interest 
concerns the stagnation points of the flow, those points of the phase space where J = 0. 


III. STAGNATION POINTS OF THE FLOW 

In m it was conjectured that the zeros of the Husimi function are also saddle points of 
J. This means that whenever a zero of the function occurs, there are topological differences 
between the classical and quantum probability flows induced by these new stagnation points, 
and the phase space dynamics of these regimes are different. In this section we prove that 
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this conjecture is true, and we analyse how the existence of the saddle points can change 
the dynamical behaviour of a system.The quantum current ([9]), like the classical one, has 
a gauge freedom under the transformation J i-> J + for any phase space valued real 

function $. Apart from this gauge, in this work we regard (§ as the probability current 
associated to the Husimi function due to its special factorization (12) shown below. 

Let Zo be a stagnation point, J (bo,~o) = 0. The probability current can be expanded 
around z 0 as 


J 


J 


dJ 

{z ~ Zo) Tz 

, i dJ 

(z “ Zo) Tz 


ZQ,Zo 


zq,zq 


t- - \ dJ 
+ (z - zo) Y _ 

, , dJ 

+ (z - z a ) m 


Z0,Z0 


20,^0 


The topological behaviour of the current about this point is determined by the eigenvalues 
A of the matrix Q of the linear coefficients of the expansion. This matrix is also the vector 
gradient of the current, and is given by 


Q [l A 


20,20 


dJ_ dJ_ 

dz dz 

dJ dJ_ 

dz dz 


20,20 


As an illustration consider the classical case, where the current is given by J c i = 
with Q — \i(}(z,z)\~. There are two sets of stagnation points, given by Q = 0 and by = 0. 
The eigenvalues of Q at these points are 


A d = 0, if Q = 0, 

A ci = ±2vaitK, if = 0; 


( 10 ) 


where K. is the Hessian matrix of the Hamiltonian function evaluated at the stagnation point: 


JC = 


d 2 H d 2 H 
dz 2 dzdz 

d 2 H d 2 H 
dzdz dz 2 


zo,zo 


When Q — 0 the structure of the stagnation point cannot be inferred. When = 0 the 
eigenvalues are either real numbers with opposite sign if the Hessian determinant is positive, 
or pure imaginary conjugate numbers if the determinant is negative. In the former case the 
point is a saddle of the flow, with an attractive and a repulsive direction; in the latter, it is 
a vortex [6]. 












For each stagnation point Zq the topological characterization of the probability current 
can be done by the index / of the flow, which counts the number of times the current rotates 
completely while we move clockwise around the point. Taking a loop L around Zq containing 
no other stagnation point, / is calculated as the clockwise integral of the angle 0 between 
the flow and some fixed reference axis in L: 

1 {Za) = ^ 

One counterclockwise rotation of J adds —1, whereas one clockwise rotation adds +1. In 
general, for saddle points I = — 1 and for vortices / = +1. In this way the index define a 
topological charge ±1 for each point (Si rrr . 

When we analyse the quantum flow, the behaviour of stagnation points is different from 
the classical flow. For pure states, the Husimi function can be factored as 


Q(z,z) = e Zz 9(z)9(z ), 


( 11 ) 


where 9 is given in Q. As the quantum current (J9]) depends only on the derivatives of Q 
with respect to z, it can be written as 

(—1 ) k d 2k+l H d l ~ l (e~ zz 9 (z)) 


^ ^ (k + l)\ dz k+l dz k 

k>0,l> 1 v ' 


dz l 


-i 


= 9 (z)f(z,z) 


( 12 ) 


Here there are two possibilities that produce J = 0: 9 = 0 or / = 0. The points that satisfy 
the condition 9 = 0, which are the zeros of the Husimi function, will be named the trivial 
stagnation points , while those that satisfy / = 0, given by an intricate relation between the 
phase space functions Q and H , will be called the non-trivial stagnation points. 

The eigenvalues of the vector gradient for both classes of stagnation points can be calcu¬ 


lated using the factorization (12): 


G = 


nd£ 9(6 f) 

u dz dz 


d ( S f) ndf 

dz U dz 

and are given by the roots of the secular equation 

A2 ~ A (4 - + i JF) +ti ¥jF-( 9 IF + 7r-') (# + f/]=•)■ 

\ oz az oz oz \ oz oz J V dz dz 


(13) 
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The solutions of (13) are different for each case considered before and are given by 


A ± = ±|f|/|, if 0 = 0; 

A 1 * 1 = | (trQ ± \/tr 2 Q — 4det Q^j , if / = 0. 


(14) 


Therefore, the trivial stagnation points are saddles of the flow, and their indices are equal 
to —1. This is the proof of the previous conjecture. There are two possibilities to the 
eigenvalues of the non-trivial stagnation points. First, when the term under the square 
root is positive, the eigenvalues are both negative (positive) numbers, and the stagnation 
point is an attractive (repulsive) node. Second, when the term inside the root is negative, 
the eigenvalues are a pair of complex conjugate numbers, and in this case the stagnation 
point is an attractive (repulsive) spiral if their real parts are negative (positive). For both 
possibilities, the real parts of the A ±; s add to — ^ and the index of the point is equal to +1. 

During the time evolution of the quantum state, the movement of the Husimi function is 
accompanied by the movement of its zeros. In view of the Poincare-Hopf theorem the total 
index of the flow must be conserved during the dynamics, thus the emergency of a saddle 
point must always be accompanied by the emergency of a non-trivial stagnation point; for 
this reason stagnation points exist in pairs. 

In HU it has been observed that the stagnation points in the pair move closely to each 
other in the phase space and form a structure similar to the one depicted in Figure 1. Since 
this structure is similar to a dipole with opposite charges, in this work we name it topological 
dipole. In the next section we analyse a toy model where the presence of the topological 
dipoles works as a signature of differences between two regimes of transmission across a 
potential barrier. 


IV. TUNNELING IN THE GAUSSIAN BARRIER 

We consider a particle of mass m scattering off a one-dimensional Gaussian barrier with 
amplitude Vo and width l/\/2 ~k. We are interested in the comparison between the classical 
and quantum transmission rates through the barrier, Tq and Tq , respectively. 

The classical Hamiltonian is 

fal - .f + Vo exp (-fa 2 ) , 


(15) 
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FIG. 1. Sketch of one topological dipole, comprised by a saddle point (left, blue spot) and a spiral 
(right, red spot). The total index around both points sums to zero. The vortex can also have an 
attractive or repulsive character. 


and the quantum Hamiltonian for the model is given by 

// = ( ( + Vfj exp {-kx 2 ), 


(16) 


In general, the classical Hamiltonian is different of the averaged Hamiltonian function defined 
in 0. and the latter contains terms of higher order in h. For the operator above, the averaged 
Hamiltonian ([6]) is given by 


H 


P TX / , 2 \ 

--b — + a l/ 0 exp (- akx ) , 

2m 4 v ' 


(U) 


1 /a 

where x and p are given by expressions (2) and a = (1 + 2 ka^) ~ is a smoothing factor. 
The classical Hamiltonian is recovered from the averaged one when h —> 0. The initial state 
was chosen as a coherent state centered at Xq and po with position and momentum widths 
cr x and cr p , respectively. The parameters were set to m = oj = 1 and h = 1/100, which 
implies a x = <r p = 1/ (l0-\/2)- We also fixed V 0 — 2, k — 3 and x 0 = —4. We also define 
Pc = \/2mVo = 2, corresponding to a classical kinetic energy equal to the barrier top. For 
further details on the numerical calculations we refer the reader to the Appendix. 

In order to quantify the relative classical-to-quantum transmission T and reflection R = 
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a 


b 




x 


FIG. 2. (a) Relative transmission and reflection coefficients for different values of p o- (b) Husirni 
function for selected times for po = 1.8. Along the arrow, time steps increase by 0.2 with initial 
time 1.7. Amplitude scale relative to the highest value, chosen to be 1. 


1 — T rates we define 


Dx — 
D r = 


T, 


c 


-Q 


T c ’ 
Rc — Rq 


R 


c 


If Dr > 0 and Dr < 0 (Dr < 0 and Dr > 0) the probability of the classical particle crossing 
the barrier is higher (smaller) than the tunneling probability of the quantum particle. There 
are two different regimes of transmission according to the average initial momentum po of the 
particle, as can be seen in Figure 2a. For p 0 ~ pc the classical transmission is greater than 
the quantum one, while if p 0 ~ pc the quantum transmission is greater. We investigated the 
structure of the quantum flow for both regimes for a state with initial momentum p 0 = 1.8 
and po = 2.1. The quantum movement of the Husimi function is shown for reference in 
Figure 2b for po = 1.8. 

Figure 3 shows the logarithmic plot of the Husimi function, highlighting the position of its 
zeros, the so called stellar representation pea mu ng. The panels show the Husimi function 
at different times when the particle is hitting the barrier. In all panels it is possible to see 
a row of zeros in front of the maximum of the Husimi function (marked by an ellipse with 
the letters RZ) in a region where classical trajectories cross the potential barrier (p > 0). 
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The presence of the zeros causes the quantum flow to be highly distorted with respect 
to the classical flow. This is seen in Figure 4, which shows the Husimi function superim¬ 
posed with the quantum current. Three zeros, which are saddles of the current, and their 
corresponding center companions, are clearly visible (marked by squares, triangles and cir¬ 
cles). The flow, that would classically go through the top of the barier to the other side, 
gets partially blocked by the topological dipoles. This dynamical feature leads to a smaller 
quantum probability of transmission compared to the classical one. 

In a similar fashion the zeros of the Husimi function and the associated stagnation points 
of the current help understand the dynamics of the transmission for po > p c , when the 
quantum transmission is larger than the classical. Figure 4b shows a few classical trajectories 
superimposed with the Husimi function for p 0 = 2.1. Once again the row of zeros is visible, 
but this time they are situated in a region near the classical separatrix. Below the separatrix, 
where the Husimi function is large and the classical trajectories are reflected back, the zeros 
distort the flow again, allowing portions of the Husimi function to cross it. Notice that 
alongside with the last zero of the row, below the separatrix, is the vortex of the topological 
pair, allowing the flow to circulate around it and move to the other side of the separatrix. 
This leads to a higher quantum than classical transmission. 

In summary, the position of the zeros relative to the Husimi function’s maximum and 
to the classical flow lines are a signature of the transmission regime. This particular model 
exhibited two particular possibilities of relative position, other systems could offer different 
situations to be explored. 


V. FINAL REMARKS 

In this work we studied the dynamics of the Husimi function and the role of its zeros 
in producing stagnation points of the corresponding phase space flow. We showed that the 
zeros are stagnation saddle points of the quantum probability flow, leading to new topological 
structures in the current when compared to the classical dynamics. These new stagnation 
points of the current are created in pairs due to topological restrictions, with indexes equal 
to +1 and —1, so that each pair behaves as a topological dipole with total index equal to 
zero. 

As an example we studied the scattering of a wave packet through a Gaussian barrier, 
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FIG. 3. Gray scale log-plot of the Husimi function for time instants T = 1.9, 2.1, 2.3 and 2.5 (see 
Figure 2b), organized in the reading direction. Black represents the absolute maximum and white 
represents zero. A row of zeros (RZ) is seen in front of the Husimi maximal values, which is framed 
by the border zeros (BZ). The outer zeros (OZ) in the external region are numerical artifacts. 


where two regimes of transmission can be identified according to the initial average momen¬ 
tum of the particle. For initial average kinetic energy below the barrier top the classical 
transmission is greater than the quantum one, whereas the quantum transmission is larger 
than the classical if the average kinetic energy is above the barrier top. In each case the rela¬ 
tive position of the topological dipoles is different, accounting for the dynamical mechanism 
behind the differences between the classical and quantum flow. When the classical trans- 
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FIG. 4. (a) Zoom in the row of zeros for T = 2.1 and po = 1-8. The black continous lines are the 
energy levels of the classical Hamiltonian, which coincides with the direction of the classical flow. 
Three topological dipoles are visible in the image, marked with square, triangle and circle, (b) 
Husimi function for po = 2.1 and T = 2.2. The black continuous lines identify the classical energy 
levels and the green curve is the separatrix. The row of zeros can be seen crossing the searatrix. 

mission is greater than the quantum, the dipoles are located in regions of the phase space 
where classical trajectories cross the barrier, partially blocking the transmission. When the 
quantum transmission is greater than the classical, the dipoles are situated near the classical 
separatrix and they offer a path for probability transmission across it, a classically forbidden 
mechanism. 

The zeros of the Husimi function have already been pointed out as signatures of quantum 
phenomena. Understanding how these zeros change the phase space probability flow adds 
new information about the dynamical mechanisms of quantum phenomena in the phase 
space. 
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APPENDIX 


Quantum dynamics of the Husimi function Given an initial coherent quantum state |^o), 
with mean position and momentum x 0 an d Po; respectively, its time evolution is obtained by 
the split-time operator method (STOM) |16| . The output of this method is the wavefunction 
'if ( x ) in the coordinate representation, and the Husimi function ([3]) is obtained through a 
convolution with the coherent state function f>* (z,z]x) = ( z\x): 


Q(z,z) = 


if (x) f * (z, z ; x ) dx 


The quantum current is then evaluated with @1: 


j = i 


d'- l Q 


ih p)?i -i 2-^i 


-1 ) k d 2k+l H 


ih ^ dz 1-1 ^ (k + /)! dz k+l dz k 
In our model the average Hamiltonian ( |T7| is not a hnite polynomial and for computational 
purpose we need to truncate the sums. If we consider an h expansion of the current, the 
contribution of order h N is obtained when l + k = N + 1. For example, the classical current 
has order h°, and we would use only the l = 1 and k = 0 term; the first quantum correction, 
of order h 1 , would use l = 2, k = 0 and l = 1, k = 1, and so on. In our calculations we went 
up to N — 10. 

The STOM routine for the wavefunction was made in the range —10.0 < x < 10.0 with 
resolution Ax = 0.0025. The time step was At = 0.01. In Figures 3 and 4 the grid of the 
Husimi function has resolution of 500 x 500 points. In the sequence of images in Figure 3 
the region where the Husimi function assumes the greatest values is surrounded by a border 
of aligned zeros, with an outlying “sea of zeros”. These zeros were robust-tested and only 
the border was observed to be accurately reproduced; the position of the zeros in the sea is 
very sensitive to the range and the resolution of the STOM. 

Classical dynamics of the Husimi function The classical evolution of the initial Husimi 


function was made integrating the trajectories generated by the classical Hamiltonian (15). 
For time instant t, the classical function is 

Qfx (t),p(t);t) = Q(x,p] 0), 

where (x(t),p(t)) is the trajectory with initial condition (x(0),p(0)) = (x,p), driven by 
Hamilton’s equations of movement: 

. dH cl . dH cl 


x = 


dp 


V = 


dx 
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In the classical case we used as initial probability density the Husimi function for this state. 
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